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Abstract 

Adiabatic quantum algorithms are often most easily formulated using many-body interactions. How- 
ever, experimentally available interactions are generally two-body. In 2004, Kempe, Kitaev, and Regev 
introduced perturbative gadgets, by which arbitrary three-body effective interactions can be obtained 
using Hamiltonians consisting only of two-body interactions. These three-body effective interactions 
arise from the third order in perturbation theory. Since their introduction, perturbative gadgets have 
become a standard tool in the theory of quantum computation. Here we construct generalized gadgets so 
that one can directly obtain arbitrary fe-body effective interactions from two-body Hamiltonians. These 
effective interactions arise from the k th order in perturbation theory. 

1 Perturbative Gadgets 

Perturbative gadgets were introduced to construct a two-local Hamiltonian whose low energy effective Hamil- 
tonian corresponds to a desired three-local Hamiltonian. They were originally developed by Kempe, Kitaev, 
and Regev in 2004 to prove the QMA-completeness of the 2-local Hamiltonian problem and to simulate 3-local 
adiabatic quantum computation using 2-local adiabatic quantum computation]!]. Perturbative gadgets have 
subsequently been used to simulate spatially nonlocal Hamiltonians using spatially local Hamiltonians 6J, 
and to find a minimal set of set of interactions for universal adiabatic quantum computationpQ. It was also 
pointed out in [5] that perturbative gadgets can be used recursively to obtain fc-local effective interactions 
using a 2-local Hamiltonian. Here we generalize perturbative gadgets to directly obtain arbitrary fc-local 
effective interactions by a single application of fc th order perturbation theory. Our formulation is based on 
a perturbation expansion due to Bloch[2]. 

A A;-local operator is one consisting of interactions between at most k qubits. A general fc-local Hamil- 
tonian on n qubits can always be expressed as a sum of r terms, 

r 

i? comp = J2 c * H * 

8=1 

with coefficients c s , where each term H s is a fc-fold tensor product of Pauli operators. That is, H s couples 
some set of k qubits according to 

Hs = Cs,l 0"s,2 ... £T s ,fc, (2) 



where each operator a s ,■ is of the form 



's,j — "SJ "S,J ) 

Mil 



(3) 



where h s j is a unit vector in K , and a s j is the vector of Pauli matrices operating on the j qubit in the 
set of k qubits acted upon by H s . 
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Figure 1: The ancilla qubits are all coupled together using ZZ couplings. This gives a unit energy penalty for each 
pair of unaligned qubits. If there are k bits, of which j are in the state |1) and the remaining k — j are in the state 
|0), then the energy penalty is j(k — j). In the example shown in this diagram, the 1 and labels indicate that the 
qubits are in the state |0001), which has energy penalty 3. 



We wish to simulate H comp using only 2-local interactions. To this end, for each term H s , we introduce 
k ancilla qubits, generalizing the technique of [3]. There are then rk ancilla qubits and n computational 
qubits, and we choose the gadget Hamiltonian as 

r r 

HSad = H * nc + A X) Vs > ( 4 ) 



where 



and 



#r= E l(i -ZsM, ( 5 ) 



^ 2 

l<i<j<k 



k 



and 



V s = 2_j c ^j <Js -J ® Xs d ( 6 ) 
i=i 

c s if j = 1 / 7 s 
1 otherwise. 

For each s there is a corresponding register of k ancilla qubits. The operators X s j and Z s j are Pauli X 
and Z operators acting on the j ancilla qubit in the ancilla register associated with s. For each ancilla 
register, the ground space of H^ nc is the span of |000 . . .) and |111 . . .). A is the small parameter in which 
the perturbative analysis is carried out. 
For each s, the operator 

Xf k = X s>1 ® X s , 2 $ . . . ® X Stk (8) 

acting on the k ancilla qubits in the register s commutes with H 6&d . Since there are r ancilla registers, 
H s&d can be block diagonalized into 2 r blocks, where each register is in either the +1 or —1 eigenspace of 
its Xf k . In this paper, we analyze only the block corresponding to the +1 eigenspace for every register. 
This +1 block of the gadget Hamiltonian is a Hermitian operator, that we label #.f d . We show that the 
effective Hamiltonian on the low energy eigenstates of _ff^_ ad approximates H comp . For many purposes this 
is sufficient. For example, suppose one wishes to simulate a fc-local adiabatic quantum computer using a 
2-local adiabatic quantum computer. If the initial state of the computer lies within the all +1 subspace, 
then the system will remain in this subspace throughout its evolution. To put the initial state of the system 
into the all +1 subspace, one can initialize each ancilla register to the state 

| +) = _L(|000...} + |111...», (9) 

which is the ground state of H^ uc within the +1 subspace. Given the extensive experimental literature 
on the preparation of states of the form |+), sometimes called cat states, a supply of such states seems a 
reasonable resource to assume. 
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The purpose of the perturbative gadgets is to obtain /c-local effective interactions in the low energy 
subspace. To quantify this, we use the concept of an effective Hamiltonian. We define this to be 



d 

H cS (H,d)^Y. E ^^\> ( 10 ) 

3=1 



where \ipi), ■ ■ ■ , \tpd) are the d lowest energy eigenstates of a Hamiltonian H, and E\, . . . , Ed are their energies. 
In section |3l we calculate H c s(H^ d , 2") perturbatively to fc th order in A. To do this, we write H £ad as 

H sad = H anc + XV (11) 

where 

r 

F anc = ^ H^ nc (12) 
s=l 

and 

v = y\v. (13) 



E 

s=l 



We consider H anc to be the unperturbed Hamiltonian and XV to be the the perturbation. We find that XV 
perturbs the ground space of H anc in two separate ways. The first is to shift the energy of the entire space. 
The second is to split the degeneracy of the ground space. This splitting arises at k th order in perturbation 
theory, because the lowest power of AV^ that has nonzero matrix elements within the ground space of H anc 
is the k th power. It is this splitting which allows the low energy subspace of £f| ad to mimic the spectrum of 

j^comp 

It is convenient to analyze the shift and the splitting separately. To do this, we define 

H cS (H, d, A) ee H cS {H, d) - An, (14) 

where n is the projector onto the support of H c b(H, d). Thus, H e g(H, d, A) differs from H c g(H, d) only by 
an energy shift of magnitude A. The eigenstates of H c g (H, d, A) are identical to the eigenstates of H c g (H, d) , 
as are all the gaps between eigenenergies. The rest of this paper is devoted to showing that, for any fc-local 
Hamiltonian H comp acting on n qubits, there exists some function /(A) such that 

H eS (Hf d , 2™, /(A)) = - fe (-^) gcomp p+ + 0(A fcU) (15) 

(fe-1)! 

for sufficiently small A. Here P+ is an operator acting on the ancilla registers, projecting each one into the 
state |+). To obtain equation [T5l we use a formulation of degenerate perturbation theory due to Bloch[2l[5], 
which we describe in the next section. 



2 Perturbation Theory 

Suppose we have a Hamiltonian of the form 

H = + XV, (16) 

where has a d-dimensional degenerate ground space £^ of energy zero. As discussed in [HEj, the 
effective Hamiltonian for the d lowest eigenstates of H can be obtained directly as a perturbation series in 
V. However, for our purposes it is more convenient to use an indirect method due to Bloch (2J E] , which we 
now describe. As shown in appendix [Bj the perturbative expansions converge provided that 

l|AV|| < J, (17) 
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where 7 is the energy gap between the eigenspace in question and the next nearest eigenspace, and || • | 
denotes the operator norirQ. 

Let \ipi), ■ ■ ■ , \^d) be the d lowest energy eigenstates of H , and let Ei, . . . , Ed be their energies. For small 
perturbations, these states lie primarily within £(°K Let 

l«j) = %)> (18) 

where Pq is the projector onto For A satisfying [T71 the vectors |ai), . . . , |ay) are linearly independent, 
and there exists a linear operator U such that 

\^)=U\a 3 ) for j -1,2,..., d (19) 

and 

U\4>) = for \(f>) e £ {Q)± . (20) 
Note that U is in general nonunitary. Let 

A = XP VU. (21) 

As shown in H] and recounted in appendix [3] the right eigenvectors of A are |ai), . . . , \ad), and the 
corresponding eigenvalues are E\, . . . ,Ej. Thus A encodes all information about the low lying eigenstates 
of the perturbed Hamiltonian. 

A and U have the following perturbative expansions. Let S l be the operator 

f E— %r in>0 

S l = I &° 1 J ' (22) 
[ -P if I = 

where Pj is the projector onto the eigenspace of with energy Ef ] . (Recall that E^ = 0.) Then 

OO 

A=Y, A(m \ ( 23 ) 

m— 1 

where 

^(m) = A m P a VS h VS h . . . VS^-'VPo, (24) 

(m-l) 

and the sum is over all nonnegative integers l\ . . . l m -i satisfying 

Zi + ... + Z m _i = m-l (25) 
Zi + ... + Z p > p 0= 1,2,..., to -2). (26) 
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Similarly, IA has the expansion 

U = P a +Y,U [m \ (27) 

m— 1 

where 

= A m ^ S h VS h V . . . VS lm VP Q , (28) 

(m) 



1 For any linear operator M, 

IMI = max |(*/>|M|tM|. 
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and the sum is over 



h + ... + l m = m (29) 
h + ... + l p > p (p=l,2,...,m-l). (30) 

In appendix lAl we derive the expansions for IA and A, and in appendix [Bl we prove that condition 1 1 71 suffices 
to ensure convergence. The advantage of the method of over the direct approach of [3] is that A is an 
operator whose support is strictly within £(°\ which makes some of the calculations more convenient. 

3 Analysis of the Gadget Hamiltonian 

Before analyzing H 6ad for a general fc-local Hamiltonian, we first consider the case where H comp has one 
coefficient c s = 1 and all the rest equal to zero. That is, 

H a * m » = <r 1 a 3 ...a k , (31) 

where for each j, aj = hj ■ &j for some unit vector fij in R 3 . The corresponding gadget Hamiltonian is thus 

iJ gad = H anc + XV, (32) 

where 

H aac = J2 \{I-Z l Z j ), (33) 

l<i<j<k 

and 

k 

V = ^ 1 ®X 1 . (34) 
j=i 

Here aj acts on the j th computational qubit, and Xj and Zj are the Pauli X and Z operators acting on the 
j th ancilla qubit. We use k th order perturbation theory to show that H cS (iJ| ad , 2 fe , A) approximates H comp 
for appropriate A. 

We start by calculating A for Hf d . For # anc , the energy gap is 7 = k — 1, and ||F|| = fe, so by condition 
1171 we can use perturbation theory provided A satisfies 

A<ij^. (35) 

Because all terms in A are sandwiched by Pq operators, the nonzero terms in A are ones in which the m 
powers of V take a state in £^ and return it to £(°\ Because we are working in the +1 eigenspace of X® k , 
an examination of equation [33l shows that £^ is the span of the states in which the ancilla qubits are in 
the state |+). Thus, Pq = I (g) P + , where P + acts only on the ancilla qubits, projecting them onto the state 
|+). Each term in V flips one ancilla qubit. To return to £ l -°\ the powers of V must either flip some ancilla 
qubits and then flip them back, or they must flip all of them. The latter process occurs at k th order and 
gives rise to a term that mimics if com P. The former process occurs at many orders, but at orders k and 
lower gives rise only to terms proportional to Pq. 

As an example, let's examine A up to second order for k > 2. 

A (<2) = X p VP + X 2 P VS 1 VP Q (36) 

The term Po^-fo is zero, because V kicks the state out of £(°\ By equation IMl we see that applying V to a 
state in the ground space yields a state in the energy k — 1 eigenspace. Substituting this denominator into 
S 1 yields 

A^ =—^PoV 2 P . (37) 
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Because V is a sum, V 2 consists of the squares of individual terms of V and cross terms. The cross terms flip 
two ancilla qubits, and thus do not return the state to the ground space. The squares of individual terms 
are proportional to the identity, thus 

A® = A 2 a 2 P (38) 

for some A-independent constant a 2 . Similarly, at any order m < fc, the only terms in V m which project 
back to are those arising from squares of individual terms, which are proportional to the identity. Thus, 
up to order k — 1, 



^-D= 5> ro A™ P 



(39) 



where the sum is over even m between zero and k — 1 and cto, a%, . . . are the corresponding coefficients. 

At fc th order there arises another type of term. In V there are fc-fold cross terms in which each of the 
terms in V appears once. For example, there is the term 



X k P (a 1 ® X 1 )S 1 {a 2 <g> X 2 )S 1 . . . S\a k ® X k )P Q 



The product of the energy denominators occurring in the S* 1 operators is 



1 



(-1)*- 



n ~j(k-j) a*- 1)!)2- 



Thus, this term is 



which can be rewritten as 



(_l)fc-iA* 
((fc-l)O 2 



P (o-i ® Xl)((72 ® X 2 ) . . . (CTfc (g) Xfe)P , 



-(-A)* 



r Po(a 1 a 2 ... ( i fc ®X® fe )Po. 



(40) 



(41) 



(42) 



(43) 



((fc-1)!) 2 

This term mimics H comp . The fact that all the S operators in this term are S 1 is a general feature. Any 
term in where h . . . h-i are not all equal to 1 either vanishes or is proportional to p). This is because 
such terms contain Pq operators separated by fewer than k powers of V, and thus the same arguments used 
for m < k apply. 

There are a total of fc! terms of the type shown in expression 1401 Thus, up to k th order 



which can be written as 



A (<k) = /(A)Po + K A) _ X® k )P Ql 

(fc - 1)! 



A (<k) = /(A )p + ^ y p (g comp g, X ® k )P 

(fc - 1)! 



(44) 



(45) 



where /(A) is some polynomial in A. Note that, up to fc th order, A turns out to be Hermitian. Up to error 
terms of order higher than fc, the effective Hamiltonian is therefore UALfi . Thus by equation l45l 



H eS {Hf d ,2 k ) = Uf(\)P U^+U 



~fc(-A) A 
(A -1)1 



P (P comp $ X m )P + 0(A fc+i ) 



w 



= /(A)n+w 

since UP W = n. Thus, 

P off (Pf d ,2 fc ,/(A))=W 



-fc(-A) fc 
(fc-1)! 

-fc(-A)* 



Po(P c 



X®*)P + 0{\ k+1 ) 



(fc-1)! 



Po(H 



comp 



.X^)Po + 0(X k+1 ) 



W 



(46) 



(47) 
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To order A fc , we can approximate U as Pq since the higher order corrections to U give rise to terms of order 
\ k+1 and higher in the expression for H eS (Hf d , 2 k , /(A)). Thus, 



P cff (Pf d , 2 fe , /(A)) = * ( _ ^ P (g comp ® A 8fc )P + 0(A fe+1 ). (48) 
Using Po = <8> P+ we rewrite this as 

P cff (Pf d , 2 fc , /(A)) = =Mz^! g com P g, p+ + 0(A fc+i). (49) 

Now let's return to the general case where H comp is a linear combination of fc-local terms with arbitrary 
coefficients c s , as described in equation [TJ Now that we have gadgets to obtain fc-local effective interactions, 
it is tempting to eliminate one fc-local interaction at a time, by introducing corresponding gadgets one by 
one. However, this approach does not lend itself to simple analysis by degenerate perturbation theory. This 
is because the different fc-local terms in general act on overlapping sets of qubits. Hence, we instead consider 



^gad = ^2 V s (50) 

8=1 

as a single perturbation, and work out the effective Hamiltonian in powers of this operator. The unperturbed 
part of the total gadget Hamiltonian is thus 

r 

P anc = HT 1C , (51) 

s=l 

which has energy gap 7 = k — 1 . The full Hamiltonian is 

^gad = H a,nc + X ygad > ^ 

so the perturbation series is guaranteed to converge under the condition 

A < * '~\ H (53) 

As mentioned previously, we will work only within the simultaneous +1 eigenspace of the X® k operators 
acting on each of the ancilla registers. In this subspace, iJ anc has degeneracy 2™ which gets split by the 
perturbation AV so that it mimics the spectrum of H comp . 

Each V s term couples to a different ancilla register. Hence, any cross term between different V s terms 
flips some ancilla qubits in one register and some ancilla qubits in another. Thus, at fc th order, non-identity 
cross terms between different s cannot flip all fc ancilla qubits in any given ancilla register, and they are 
thus projected away by the Po operators appearing in the formula for A. Hence the perturbative analysis 
proceeds just as it did when there was only a single nonzero c s , and one finds, 

P cff (Hf d , 2", /(A)) = ^Ew Po c ° H ° X ® k ) Po + °( xk+1 )> (M) 



\s=l 



where Xf k is the operator X® k acting on the register of fc ancilla qubits corresponding to a given s, and 
/(A) is some polynomial in A of degree at most fc. Note that coefficients in the polynomial /(A) depend on 
P comp . As before, this can be rewritten as 



H cS (Pf d , 2", /(A)) = ^L^! g comp (g>p+ + 0(A fc+i )j , - - , 
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where P + acts only on the ancilla registers, projecting them all into the |+) state. Hence, as asserted 
in section [T] the 2-local gadget Hamiltonian H 6ad generates effective interactions which mimic the A:-local 
Hamiltonian # com P. 

For a polynomial time adiabatic quantum computation one needs a Hamiltonian that varies smoothly in 
time and has an eigenvalue gap at worst polynomially small. Let H (t) be a fc-local Hamiltonian of this type. 
For each time t one can construct the corresponding instantaneous gadget Hamiltonian H s&d (t) as described 
in section[T]and equation[53l It is not hard to show that H gad (t) varies smoothly in time and has a gap that is 
polynomial in n for any fixed k. Thus H md (t) is a 2-local polynomial-time adiabatic algorithm that simulates 
the original /c-local algorithm H(t). In addition to adiabatic quantum computation we expect that fc th order 
gadgets may have many other applications in quantum computation, such as proving QMA-completeness. 



4 Numerical Examples 

In this section we numerically examine the performance of perturbative gadgets in some small examples. As 
shown in section [31 the shifted effective Hamiltonian is that given in equation [551 We define 



H 



id 



-fc(-A) fc 

(k-iy. 



H comp ~ p 



(56) 



.ffeff consists of the ideal piece H ld , which is of order X k , plus an error term of order A fc+1 and higher. For 
sufficiently small A these error terms are therefore small compared to the H ld term which simulates iJ com P. 
Indeed, by a calculation very similar to that which appears in appendix [Bl one can easily place an upper 
bound on the norm of the error terms. However, in practice the actual size of the error terms may be 
smaller than this bound. To examine the error magnitude in practice, we plot - — || gid |° m figure [2] using 
direct numerical computation of H c s without perturbation theory. /(A) was calculated analytically for these 
examples. In all cases the ratio of ||-ff ld — H e g || to ||i? ld || scales approximately linearly with A, as one expects 
since the error terms are of order X k+1 and higher, whereas H ld is of order A fe . 

0.25 . 1 , , , 



llff* 1 - -ge ff || 



0.2- 



0.15 



o 


XYZ+XYY 


□ 


XYZZ 


A 


XYZ 



0.1 



0.05- 



0.05 



0.1 



A 



-k(-X) k 



(fc-1)! 



H colnp is plotted. We examine 



Figure 2: Here the ratio of the error terms to the ideal Hamiltonian H' d 

three examples, a third order gadget simulating a single XYZ interaction, a third order gadget simulating a pair 
of interactions XYZ + XYY, and a fourth order gadget simulating a fourth order interaction XYZZ. Here H e s is 
calculated by direct numerical computation without using perturbation theory. As expected the ratio of the norm of 
the error terms to H' d goes linearly to zero with shrinking A. 
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A Derivation of Perturbative Formulas 

In this appendix we give a self-contained presentation of the derivations for the method of degenerate 
perturbation theory used in this paper. We closely follow Bloch[2]. Given a Hamiltonian of the form 

H = # (0) + XV (57) 

we wish to find the effective Hamiltonian induced by the perturbation XV on the ground space of H W . In 
what follows, we assume that the ground space of H 1 - ' has energy zero. This simplifies notation, and the 
generalization to nonzero ground energy is straightforward. To further simplify notation we define 

V = XV (58) 

Suppose the ground space of is d-dimensional and denote it by £(°\ Let \ipi), ■ ■ ■ , \ipd) be the 
perturbed eigenstates arising from the splitting of this degenerate ground space, and let E±, . . . , Ed be their 
energies. Furthermore, let |ay) = Po\ipj) where Pq is the projector onto the unperturbed ground space of 
H^K If A is sufficiently small, . . . , \ad) are linearly independent, and we can define an operator IA such 
that 

and 

Now let A be the operator 



U\a s ) = 1^-) (59) 
U\cj>) = M\<t>) e £ (0)± . (60) 
A = P VU. (61) 



A has \ai), . . . , \ad) as right eigenvectors, and E\, . . . , Ed as corresponding eigenvalues. To see this, note 
that since has zero ground state energy 

P V = P {H^+V)=P H. (62) 

Thus, 

A\ccj) = P VU\ aj ) 

= Pt>V\4>j) 

= iWi) 

= P6Ej\4>j) 

= Ejlatj). (63) 

The essential task in this formulation of degenerate perturbation theory is to find a perturbative expansion 
for IA. From IA one can obtain A by equation 1611 Given A, one can easily claulate its right eigenvectors 
\ati), ■ ■ ■ , \ctd) an d the corresponding eigenvalues E±, . . . , Ed- Then, by applying U to \aj) one obtains 
So, given a perturbative formula for hi, all quantities of interest can be calculated. 
To derive a perturbative formula for IA, we start with Schrodinger's equation: 

H\iP j ) = E j \i> j ). (64) 
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By equation [511 left-multiplying this by Pq yields 

P Q V\^ j )=E j \a j ). (65) 

By equation 1501 

UP = U. (66) 

Thus left-multiplying equation [55] by U yields 

UV\^)=Ej\^j). (67) 

By subtracting [571 from [Ml we obtain 

{H-UV)\ipj) =0. (68) 
The span of \tpj) we call £ . For any state \/3) in £ we have 

(H-UV)\/3) = 0. (69) 

Since £ £ for any state I7), it follows that 

(H - UV)U = 0. (70) 

This equation can be rewritten as 

H {0) U = -VIA + UVU. (71) 

Defining Qq = 1 — Pq we have 

U = P Q U + Q U. (72) 
Substituting this into the left side of [7T1 yields 

H m Q Q U = -VU+UVU, (73) 

because H^P = 0. In £^ ± , has a well defined inverse and one can write 

QoU = -j^Q (VU - UVU). (74) 

Using equation [721 one obtains 

IA = P U- j^Qa(VU - UVU). (75) 
By the definition of U it is apparent that PqU — Pq, thus this equation simplifies to 

U = P Q -j^ Ql Qo(VU-UVU). (76) 

We now expand U in powers of A (equivalently, in powers of V), and denote the m th order term by U^ m >. 
Substituting this expansion into equation [76] and equating terms at each order yields the following recurrence 
relations. 



U {0) = P (77) 

rn— 1 



u(m) = -lk) Qo 



P =i 

Note that the sum over p starts at p = 1, not p — 0. This is because 



(m = 1,2,3...) (78) 



^ m QoU^ = j^ J Q P = 0. (79) 
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Figure 3: From a given m-tuple (hjfo, ■ ■ ■ ,lm) we construct a corresponding stairstep diagram by making the j* 
step have height lj, as illustrated above. 



Let 

S l = { (-^ (0) ) . (80) 

[ -P ifZ = 

£/( m ) j s f the form 

U (m) =J2's h VS h V...S lm VP , (81) 
where ^2 is a sum over some subset of m-tuples (h 7 h, ■ ■ ■ Jm) such that 

^ > (t = l J 2,...,m) (82) 

Zi + Z 2 + • • • + lm = m. (83) 

The proof is an easy induction. clearly satisfies this, and we can see that if has these properties 
for all j < to, then by recurrence [78j also has these properties. 

All that remains is to prove that the subset of allowed TO-tuples appearing in the sum ^ are exactly 
those which satisfy 

h + ... + l p >p {p= 1,2,..., m- 1). (84) 

Following [2], we do this by introducing stairstep diagrams to represent the TO-tuples, as shown in figure[3] 
The m-tuples with property [84] correspond to diagrams in which the steps lie above the diagonal. Following 
[2] we call these convex diagrams. Thus our task is to prove that the sum is over all and only the convex 
diagrams. To do this, we consider the ways in which convex diagrams of order to can be constructed from 
convex diagrams of lower order. We then relate this to the way is obtained from lower order terms in 
the recurrence [751 

In any convex diagram, 1\>1. We now consider the two cases l\ = 1 and > 1. In the case that l\ = 1, 
the diagram is as shown on the left in figure 2] In any convex diagram of order to with li — 1, there is an 
intersection with the diagonal after one step, at the point that we have labelled c. The diagram from c to b 
is a convex diagram of order to — 1 . Conversely, given any convex diagram of order to — 1 we can construct 
a convex diagram of order to by adding one step to the beginning. Thus the convex diagrams of order to 
with l\ = 1 correspond bijectively to the convex diagrams of order to — 1. 

The case l\ > 1 is shown in figure [4] on the right. Here we introduce the line from a' to b' , which is 
parallel to the diagonal, but higher by one step. Since the diagram must end at b, it must cross back under 
a'b' at some point. We'll label the first point at which it does so as c'. In general, c' can equal b' . The curve 
going from a' to c' is a convex diagram of order p with 1 < p < to — 1, and the curve going from c to b 
is a convex diagram of order n — p — 1 (which may be order zero if c! = b'). Since c' exists and is unique, 
this establishes a bijection between the convex diagrams of order to with l\ > 1, and the set of the pairs of 
convex diagrams of orders p and n — p — 1, for 1 < p < n — 1. 

Examining the recurrence 1781 we see that the l\ = 1 diagrams are exactly those which arise from the 
term 

QO y U ( m -l) (85) 
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h > 1 



b' b 




Figure 4: A convex diagram must have either h = 1 or h > 1. In either case, the diagram can be decomposed as a 
concatenation of lower order convex diagrams. 



and the l± > 1 diagrams are exactly those which arise from the term 



Qo_ 



(86) 



0=1 



which completes the proof that J2' is over the m-tuples satisfying equation [84] 

B Convergence of Perturbation Series 

Here we show that the perturbative expansion for U given in equation [27] converges for 

w\\<\- 



(87) 



By equation [22 the convergence of U also implies the convergence of A. Applying the triangle inequality to 
equation [37] yields 



(m)i 



||W||<l + EllW 



Substituting in equation [281 and applying the triangle inequality again yields 

oc 

\\U\\ < 1+ £ A m ^||5 (l ...VS'-Woll- 

m=l (m) 

By the submultiplicative property of the operator norm, 



(89) 



iiwii < i + e xm E ii^ii • wn ■ ■■ wn ■ \\s lm \\ ■ wn ■ ii^i 

m=l (m) 



| Pol I = 1, and by equation |2"21 we have 



Since the sum in equation 1901 is over l\ 



\S l \ 



1 1 



( £ (o)y y 



Z rn = m, we have 



iiwii <i+ EE 



\\\Vf 



(90) 



(91) 



(92) 



m=l (rn) 
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The sum J2( m ) 1S over a subset of the m-tuples adding up to to. Thus, the number of terms in this sum is less 
than the number of ways of obtaining m as a sum of m nonnegativc integers. By elementary combinatorics, 
the number of ways to obtain n as a sum of r nonnegative integers is f n+r thus 

„ „ ^ (2m - 1\ \\\V\\ m , s 



Since 



we have 



2 g 3 /2TO-l\ =22m _ 1; (94) 



3=0 



J 



'2m- 1 

TO 

Substituting this into equation 1931 converts it into a convenient geometric series: 



< 2 2m -\ (95) 



W|| <l+^2 2 '»- 1 ^p. (96) 

m— 1 



This series converges for 
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